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The magnetic interactions of iron superconductors are studied within the Hartree-Fock (HF) and 
Heisenberg descriptions of the same five orbital model. The experimentally observed (tt, 0) order 
is more stable in the HF approximation than in the Heisenberg one. Due to virtual transitions 
involving filled orbitals, the first and second nearest neighbor exchange constants Ji and J2 of the 
strong-coupling description become ferromagnetic (FM) at large Hund's coupling. Ferromagnetic 
tendencies are enhanced with electron doping and suppressed with hole doping, the electron-hole 
asymmetry being less important for small interactions. The Hartree-Fock phase diagram lacks a 
ferromagnetic phase but includes a double stripe (DS) state. In undoped compounds this state 
is unstable towards phase separation (PS), but it can be stabilized with electron-doping. In both 
descriptions, (tt, tt) antiferromagnetic order appears on the phase diagram correlated with an unusual 
orbital reorganization. A non-monotonous dependence of Ji and J2 on the Fe-As-Fe angle is found. 
The implications of our findings for the magnetic state of iron superconductors are discussed. 

PACS numbers: 75.10.Jm, 75.10.Lp, 75.30.Ds 
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The origin of the magnetism in iron superconductors, 
believed to be key to explain their high-Tc supercon- 
ductivity, is still unsettled. Most iron pnictides order 
with (tt, 0) momentum, AF in the a;-direction and FM 
in the y-direction[l|,|3|. The itinerant versus strong cou- 
pling origin of magnetism is at the heart of the debate. 
The itinerant picture relies on the metallicity of the AF 
and an approximate nesting between hole and electron 
pockets [3|-l5(. In the strong coupling description, local- 
ized spins interact AF with their first and second near- 
est neighbors [6|,|7|. Clasically, a (tt, 0) state arises when 
J2 > Ji/2. A large J2 was justified by the As-mediated 
exchange between Fe atoms. In contrast, FeTe orders 
FM along one of the diagonals and AF along the other 
one in a double stripe pattern. A large exchange interac- 
tion to third nearest neighbors J3 has to be introduced to 
explain the stability of this state within a localized pic- 
ture [Sj . A relation between its magnitude and the Fe-Te 
layer geometry has been claimed [9| . No nesting features 
compatible with this ordering have been observed [lO|. 

Little is known about the value of the exchange con- 
stants. In strongly correlated multiorbital systems, the 
interplay between the orbital filling and magnetism gives 
rise to a variety of magnetic phases. In iron superconduc- 
tors, estimates for Ji and J2 have been restricted to a few 
ab-initio calculations for specific compounds [8, 9.,lll-ll3|. 
The extraction of the exchange parameters from neutron 
experiments is still controversial 14l4l6l|: An unexpected 
FM Ji has been recently proposed to fit the spin-wave 
spectrum of iron chalcogenides [16| . 

In this letter we analyze the magnetic interactions 
of iron superconductors on the basis of a five orbital 
model 17| treated both within HF and Heisenberg ap- 
proximations. We focus on the AF or FM character of 
the strong coupling exchange constants Ji and J2 , on how 
they depend on the doping and on the Fe-As layer geom- 
etry, and check whether (tt, 0) is the actual ground state. 



We find that in undoped compounds Ji and J2 have a 
non-trivial dependence on the Fe-As-Fe angle and may 
become FM at large Hund's coupling Jh due to virtual 
transitions involving filled 3z^ — r^ or x'^ — y^ orbitals. 
FM tendencies are suppressed with hole-doping, but en- 
hanced with electron-doping. The region of stability of 
the (tt, 0) state increases in HF with respect to the strong 
coupling predictions. An unusual orbital reorganization 
stabilizes (tt, tt) ordering in both HF and Heisenberg de- 
scriptions for a wide range of parameters. A DS state 
shows up at large U and Jh in the HF phase diagram. 
This DS state is charge modulated and unstable towards 
phase separation for n = 6 but it can be stabilized with 
electron doping. 

The 5-orbital Hamiltonian includes intraorbital U and 
interorbital [/', Hund's Jr and pair hopping J' terms. 
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i,j label the Fe sites in the Fe unit cell, a the spin and 



7 and /3 the five Fe d-orbitals yz, zx, xy, 3z 



and 



X — y , with X and y axis along the Fe-Fe bonds. The 
hopping amplitudes, restricted to first and second neigh- 
bors, depend on the angle a formed by the Fe-As bonds 
and the Fe-plane[l7|. a = 35.3°, corresponding to the 
regular Fe-As tetrahedra, is assumed except when in- 
dicated^ We use [/' = [/- 2Jh [11, J' = Jh, and 
take 



17| 



ckj/ 



= 0.02, e 



zx,yz 



= 0, 



£32 



— —0.55 and 



ea.2_y2 — —0.6 for the crystal field. Energies are in units 
of (pda)"^ /\ed — Cpl '^ 1 eV, with pda the a overlap be- 
tween the Fe-d and As-p orbitals and |erf — ep| their energy 



difference [17[. For details on the mapping to a classical 
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FIG. 1: (Color online) (a) Sketch of the two different orbital configurations {X = x ~ y 
for n = 6. (b) Exchange constants for n = 6, U — 5 eV, and a = 35.3° as a function of Hund's coupling, (c) Illustration of the 
virtual transitions involving filled orbitals responsible for the FM exchanges for n = 6 at large Jh/U. (d) Exchange constants 
as a function of Jh/U for n — 5 and U — 5 eV assuming a S = 5/2 state. A (tt, tt) ground state is always favored, (e) Energies 
of the magnetic ground states as a function of Jh/U for the same parameters as in (b). The double stripe order, although 
never the ground state, is included for comparison. Vertical lines in (b) and (e) separate regions with different ground states 
[(7r,7r),(7r,0) and FM]. 



Heisenberg model and HF treatment see the supplemen- 
tary information (SI) [l9l. 

To second order in perturbation theory, starting from 
localized atomic states, the 5-orbital Hamiltonian is 
mapped onto a classical Heisenberg Hamiltonian: 
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and 



with Si — J^B^i^fi t^'2 atomic moment and («,j) 
((i,j)) restricted to first and second nearest neighbors 
respectively. We focus on the largest possible spin state, 
with each orbital being half-filled or filled. For n = 6, as 
in undoped compounds, this corresponds to S* = 2. This 
large spin state dominates the (tt, 0) mean field phase 
diagram at large U 20|, [21 [ . Due to the small crystal field 
splitting, we consider two possible atomic configurations, 
labelled by {A}, for the S = 2 state with filled x^ — y^ or 
3^2 _ j,2 orbitals, see Fig. [T] (a). 

As expected, the exchange constants in both atomic 
configurations decrease monotonically with Hund's cou- 
pling, see Fig. [Ijb). But the decrease is notably steep, 
mainly for Ji and at large J^/U , where a change in the 
slope happens and the exchange constants even become 
FM. This strong dependence at large Jh is due to vir- 
tual transitions from a filled orbital to a half-filled one 
on a neighbor atom, see Fig. [TJc) and SI [l^. These 
transitions are favored by the small crystal field splitting 
characteristic of iron pnictides. On the other hand, there 
are no filled orbitals at n = 5 in its highest spin state 
{S = 5/2), implying a weaker dependence on Jh/U of 
the exchange parameters, always AF, see Fig.[TJd). The 
opposite trend is found for n = 7, not shown, with two 
completely filled orbitals and enhanced FM behavior. 

For n = 5, J2 < Ji/2 and (tt, tt) ordering is found. For 
n = 6, in the x'^ — y^ configuration, favored by the crystal 



field splitting, the (tt, O) state is lowest in energy, except 
at large Jh/U for which FM is expected. When com- 
paring the energies between the ordered states in both 
orbital configurations in Fig. [TJe), the (7r,7r) order with 



filled 32 



becomes the ground state in a wide range 



of parameters. The selection of (tt, tt) in the 3z — r con- 
figuration is in accordance with the exchange constants 



t3z' 



tSz-' 



ratio Jj^ ~^ < Ji^' ~^ /2. Remarkably, the gain in mag- 
netic energy, including Eq, compensates for the cost in 
crystal field {e^^^-r^ — ex^-y^ = 0.05), see Sl[l9|. 

The appearance of (tt, tt) as the ground state is repro- 
duced in the HF phase diagram in Fig. [2] As in the 
strong coupling picture, (tt, tt) in a 3z^ — r^ configuration 
and (tt, 0) in a x^ — y^ configuration dominate the phase 



diagram. In the (tt, 0) state, 3z 



is emptied while 



X - 
3z2 



y gets filled as U increases [20|- In the (tt, tt) state, 
— r^ fills with increasing U while the other orbitals 
tend to half filling, see SI [l9|. The transition between 
(tt, tt) and (tt, 0) as Jh increases is accompanied by a 
charge transfer between x'^ — y^ and 3z^ — r^ . 

The (tt, 0) state is more prominent in the HF phase di- 
agram than in the Heisenberg description, specially close 
to the NM region. As previously found [20j , two different 
(tt, 0) states show up at the HF level. In the LM state, 
opposite orbital magnetizations within the same atom 
result in a low magnetic moment which violates Hund's 
rule. This state is stabilized thanks to the anisotropy 
of the orbital exchange constants [20| and has been pro- 
posed to explain the low magnetic moment found exper- 
imentally 0, m [ii]. In the HM state, all the orbital 
magnetizations point in the same direction. The strong 
coupling predictions, dashed lines in Fig. [31 are valid for 
comparison with this HM state. 

At large Jh the kinetic energy effects present at the 
HF level prevent the FM solution found in the strong 
coupling description to arise. Instead, a double stripe 




U(eV) 
FIG. 2: (Color online) Hartree Fock magnetic phase diagram 
calculated at fixed density n = 6 as a function of U and 
Jh /U . Grey stands for non-magnetic (NM), blue for the AF 
(tt, 0) high moment (HM) state satisfying Hund's rule, red 
for a the AF (tt, 0) low moment (LM) state which violates 
Hund's rule, orange for AF (tt, tt) and green for the double 
stripe (DS) state with charge modulation (see text). Shaded 
areas mark insulating phases. Superposed dashed lines on the 
magnetic regions show the phase transition lines between the 
(tt, tt) and HM (tt, 0) states and between the HM (vr, 0) and 
the FM phase predicted by the S = 2 Heisenberg model. 
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FIG. 3: (Color online) Mean field magnetic phase diagram 
in the grand canonical ensemble as a function of U and Jh /U 
and for n = & {&) and as a function of U and n for Jh /U — 
0.22 (b). Same color code as in Fig. ([2]) applies. White regions 
are unstable towards phase separation. 



solution is found at large U and Jh- This DS state is 
charge modulated (5 or 7 electrons for the n — 6 case) 
in a checkerboard fashion. Note that in the Heisenberg 
picture, a homogeneous DS state, though not the ground 
state, is lower in energy than both (7r,7r) and (7r,0) at 
large Jh/U, see Fig. [TJ^e). The homogeneous DS state, 
with n electrons in each atom, is difficult to stabilize in 
the HP calculations, mainly at large Jh- When stabi- 
lized, it has larger energy than the charge modulated DS 
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FIG. 4: (Color online) (a) First and second nearest neighbor 
exchange constants Ji and Ji as a function of a, the angle 
formed by the Fe-As bond and the Fe plane, for n = 6, U — 5 
eV, and Jh/U = 0.22. (b) Energies of the most competitive 
magnetic orderings as a function of a for the parameters in 
(a). The sequence of magnetic orderings with increasing a 
is (tTjO), (tTjTt), and FM. Vertical lines delimitate the experi- 
mentally relevant values of a. 



state. As shown in Fig. [31 the DS state is unstable to- 
wards phase separation at n = 6, see SI [19|, but can 
be stabilized by electron doping. As expected for first 
order transitions, PS appears at the boundary between 
the different magnetic phases. The instability at n > 6.5, 
which only involves the DS state, has a different origin: 
it is caused by a negative compressiblity related to the 
charge modulation. 

Fig. ^h) also evidences an electron-hole asymmetry 
around the undoped composition n = 6. In agreement 
with the strong coupling predictions, for n = 5 (tt, tt) 
correlations dominate, even at large Jh (calculations at 
Jh/U = 0.28 give similar results). When n increases, the 
(tt, O) state becomes lower in energy, being the ground 
state in a wide region of parameters. For n > 6 and in- 
termediate or large U the DS state is found. A larger 
tendency to DS solutions with increasing n is also ex- 
pected in the Heisenberg description. Within HF, the 
electron-hole asymmetry close to n = 6 gets reduced as 
U decreases. 

We finally analyze the effect of the Fe-As-Fe geome- 
try on the magnetic interactions. A direct relation be- 
tween the Fe-As-Fe angle and the critical superconduct- 



mg temperature 



25| as well as on the magnetic 



ordering [16[ has been claimed. As illustrated in Fig. 21 
the exchange constants vary non-monotonically with a 
due to the relative importance of the hopping (to both 



first and second neighbors) via the As [17[. The depen 



dences of Ji and J2 on a are different for the x^ — y^ and 
3z^ _ j,2 configurations, with J1/J2 changing with a. For 

2_ 2 

some values of a, Ji ( Jf ^ in special) and J2 are FM. 
An alternance of magnetic states involving orbital reor- 
ganization, similar to the one discussed above, is found. 
For the value of Jh/U chosen, the transitions between 
these orders occur within the range of physical values of 
a, partly due to Eq, see SI [l9{ . 

In summary, we have analyzed the magnetic interac- 



tions and ground state of a five orbital model for iron 
pnictides by means of a strong-coupling Heisenberg and 
a Hartree Fock approach. The tendency to (tt, 0) order- 
ing is enhanced in the mean field description, both in the 
large Jh region and in the small U region. The kinetic 
energy nesting properties of the Fermi surface may be im- 
portant for the latter area, where HF is expected to give 
a better description and the resistivity shows anisotropic 
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behavior consistent with experiment [2l|. The value and 



sign of the exchange constants of the strong coupling 
model depend non-trivially on the Fe-As-Fe angle, the 
orbital filling, the number of electrons per Fe, and the 
Hund's coupling Jh- The small crystal field splitting 
characteristic of iron pnictides produces a strong decrease 
of the exchange constants with Jh and stabilizes an or- 
bitally reorganized (tt, tt) state in a wide range of parame- 
ters, both for the weak and strong coupling descriptions. 
In the Heisenberg description, virtual transitions from 
filled to half-filled orbitals, of little importance in many 
materials because of large crystal field splittings, induce 
a strong tendency to FM at large Jh ■ FM is absent at the 
HF level, while a double stripe solution is found at large 
Jh/U. The system shows a strong tendency towards 
phase separation in the mean field approach, specially at 
large U and Jh /U , including a region with phase sepa- 
ration between (tt, 0) and double stripe orderings. The 
double stripe state is charge modulated. Whether it is 
related to the state found in FeTe is at present not known. 

Further conclusions can be obtained from this anal- 
ysis. First, the exchange constants in the (tt, 0) and 
(tt, tt) states may be different. Thus estimating these con- 
stants by comparing the energy of the different magnetic 
states with those predicted by a mapping to a Heisen- 
berg model [Sl, y, ll2, llSj with the same Ji and J2 may 
lead to errors. Second, the orbital reorganization which 
accompanies the (tt, 0)-(7r,7r) transition implies an unex- 
pected sensitivity of the ground state to the crystal field 
splitting £3^2 _^2 — e^2_y2. Most of the proposed tight 
binding models have very similar values for this splitting 
and similar physics is expected, although the quantita- 
tive details may depend on the parameters. Third, our 
results indicate that the filled orbitals are not inert for 
magnetism, so the use of 4-orbital models which neglect 
them to estimate the exchange parameters may be ques- 
tionable [23. 

Our results stress the relevance of Hund's coupling, 
also emphasized by other authors ll|, 1261428 1 . At small 
values of Jh^ Hund's rule may be violated resulting in a 
(tt, 0) low magnetic moment state[20j. (tt, 0) ordering is 
also favored at large Jh- A large Jh produces a strong 
decrease of the exchange constants, more important with 
electron doping, leading to small AF or even FM Ji in the 
experimental range of parameters. These results could 
be relevant to understand the neutron scattering experi- 
ments. 
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selected, while Epj^'^p^.- ■ is the magnetic contribution. 



p{A} 



pW 



E^f' enters in the value of the constant E^ ' but cancels 
out in the expression of the exchange constants J^ ^ ■ We 



simplify e"^^ to El, 



-y =0and£;3. 



e3^2_,.2 



ea;2_j,2 per atom, as the other terms drop out from the 
calculation. For n — & 



Here we provide the details of the calculations we refer 
to in the text. Starting from localized atomic states we 
map the 5-orbital interacting Hamiltonian of Eq.(l) in a 
classical Heisenberg Hamiltonian with exchange interac- 
tions to first and second nearest neighbors 
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Here A labels the considered orbital configuration (as 
defined in Fig. 1 (a)). Si = ^aSi^p the atomic mo- 
ment and (j, j) and {{i,j)) are restricted to first and sec- 
ond nearest neighbors respectively. We restrict to atomic 
states with maximum total spin S {S = 2 for n = 6 and 
S = 5/2 for n = 5) and maximum jS'^j (Sz — ±2 and 
Sz = ±5/2 for n = 6 and n = 5 respectively). With this, 
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Eq ' and J{2 are 
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with i,j two first (second) nearest neighbors for Ji (J2). 

^FM-i j ^^'i ^AF-i j ^^6 ^^^ energies corresponding to 
two parallel or antiparallel spins at i and j. The factor 
of 4 (instead of 2) corrects the bond double-counting in 
Eq.dSTJ. When using Eqs. (lS2| and dSS]) we are neglect- 
ing quantum fluctuations. In this sense we are treating 
the spin classically. As discussed in the text and below, 
the classical approximation provides good understanding 
on the magnetic interactions dominant in iron pnictides. 
Quantum fluctuations are expected to be of little impor- 
tante for large spins. Spin states S = 2 and S — 5/2 
discussed in the text are large enough to justify the clas- 
sical treatment of spins. 

FM and AF energies can be written in terms of atomic 
and magnetic energies 
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where E^, ' accounts for the interaction and crystal field 
energy to charge the atoms to the n = 6 or n = 5 state 
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Here r] and v label half-filled orbitals and A refers to the 
filled one; if^ are the hopping amplitudes between or- 



bitals l3 and 7 in atoms i and j respectively 17| . AE± = 

e^-£x±{iex-e^)^ + J'^)'^^^. The first term in Eq. (|S5|) 
includes virtual transition from half-filled orbitals in one 
atom to half-filled orbitals in a neighbour atom and takes 
into account that the intermediate state is not an eigen- 
state of the pair hopping operator. The second term is 
associated with virtual transitions of an electron from 
the filled orbital to a half-filled orbital in the other atom. 
Only the electron whose spin is opposite to the magnetic 
moment of the neighbor atom can hop. The state left 
behind is not an eigenstate of the Hund's term. 

When the magnetic moments of the two atoms are par- 
allel the only transitions allowed involve the filled orbital. 
The contribution of these transitions is 
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For n = 5 and S = 5/2 Ep-f.i^:"""^ = and 
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In this latter case, the exchange constants are clearly 
AF. Whether a (tt, tt) or (tt, 0) state is favoured depends 
on the relative values of Ji and J2. As shown in Fig. 1, 
(tt, tt) ordering is always preferred. 

The n = 6 case is more interesting. For Jh — the 
contribution of the transitions which depart from filled 
orbitals is the same with FM or AF ordering. Conse- 
quently they do not affect the magnetic ordering, and 
the exchange constants, driven by the transitions from 
half-filled orbitals are AF. With increasing Jh their con- 
tribution is larger in the FM case. At a given Jh these 
transitions favor FM ordering. For small Jh the transi- 
tions from half-filled orbitals still dominate, but at large 
Jh Ji and J^ become ferromagnetic. Due to different 
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FIG. SI: (Color online) Dependence of El "^ and El"" -'' 
for n = 6 and U = 5 eV as a function of the Fe-As-Fe angle a 
for Jh/U = 0.22 (a) and Hund's coupling for a = 35.3°(b). 
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FIG. S2: (Color online) (a) and (b) Dependence of the orbital 
filling as a function of the intraorbital interaction for Jh/U = 
0.15 and Jh/U = 0.22 respectively. Vertical lines separate 
regions with different ground states. 



t] ^ in each of these terms the value of Jh at which this 
happens is different for each exchange constant. 



Another interesting aspect regards Eq . The contri- 
bution of Epj^'™p^.^ ■ to Eq ^ can be seen as the energy 
gain associated to creating magnetic moments, even if 
on average they do not order. When mapping to Heisen- 
berg models it is usual to disregard £'o because if a single 
atomic state is selected its value does not affect the dif- 
ferences in energy between magnetic states. However, for 



6, the small crystal field splitting between x 



and Sz^ — r^ requires the inclusion of two different atomic 
states. Due to different hopping amplitudes involving 

these orbitals, E^ "^ ^ -E'o^^"'^^ is found. This can 
be observed in Fig. EI] as a function of the Fe-As-Fe 
and Hund's coupling. Consequently Eq ^ helps select 
the atomic state {A}, with the a priori unexpected re- 
sult that the atomic state in which 3z^ — r^ is filled is 
preferred in a wide range of parameters. Note that the 
formation of magnetic moments is favoured in the exper- 
imentally relevant range of Fe-As-Fe angle and at large 
Hund's coupling. The large values of Eq also explain 
the large values of the total energy, despite the relatively 
small values of the exchange constants. 

We have compared these predictions with a self- 
consistent mean-field Hartree Fock calculation which in- 
cludes non-magnetic (NM), FM, and AF states with 
Q = (tt, 0), Q = (tt, tt), and DS ordering. In the DS cal- 
culation the system is divided into two interpenetrating 
lattices coupled via first nearest neighbor hopping terms. 
The axis and orbital basis is rotated and a (tf, 0) state 
along the Fe-diagonals is assumed. In the mean-field cal- 
culation only spin and orbital-diagonal average terms are 
kept: 
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We have checked that this approximation does not have 
any influence on the mean field results. Terms disre- 
garded are of order ~ 10"'*. Several initial configura- 
tions for n^ and m.y are probed for each ordered state. 
A self-consistent solution is considered to be ordered (in- 
sulating) when the magnetic moment m (gap) is larger 
than 0.001. Fig. 2 shows the phase diagram correspond- 
ing to n = 6 calculated at fixed density. The mean 
field ground state evolves from a non-magnetic solution 
with zero magnetic moment and orbital fillings close to 
the non-interacting ones at small U to (tt, 0) and (tt, tt) 
states with m = A^lb, with jj,b a Bohr magneton, as in 
an atomic S = 2 state at large U . As seen in Fig. [S2] 
every orbital except x^ — y^ tends to half filling at large 
[/, x^ — y^ becoming completely filled, in the (tt, 0) state. 
This orbital filling is expected on the basis of the crystal 
field splitting. On the other hand, in the (tt, tt) state, 
it is the 3z^ — r^ orbital the one which fills completely 
while x^ — y^ tends to half-filling. This orbital filling 
competes with the crystal field, in agreement with the 
strong coupling predictions, see main text. 

The phase diagram in Fig. 2 has been obtained com- 
paring the energies of the mean-field solutions at fixed 
density. Different phases have a different energy E versus 
density n relation. Jumps of the chemical potential \x ap- 
pear at first-order boundaries between ordered phases. In 
such situations phase separation is expected. In a range 
of densities the system decreases its energy by allowing 
relative fractions of the two neighboring magnetic phases, 
n = 6 may belong to this range of densities for values of 



S3 



the interactions close to phase boundaries. To determine 
the region of phase space unstable towards phase sepa- 
ration in Fig. 3 we work in the grand-canonical ensemble 
(fixed chemical potential) and compare E — /in for the 
different states. The chemical potential can be calcu- 
lated as dE/dn or by looking at the energy of the last 
state occupied in metallic systems, both methods giving 
the same value, /i jumps discontinuously when cross- 
ing a gap. According to Maxwell's construction phase 
boundaries between two phases 1 and 2 are given by 
El — fxni = E2 — fJ.n2. Neither phase 1 nor phase 2 
are stable in the range ni < n < n2- A phase separated 



mixture of both phases arises. 

Phase separation for a given phase will also appear 
with negative compressibility , i.e. /i decreases with in- 
creasing n. This is observed in the DS state in several 
ranges of n, including n = 6 and n < 6.5. In this case 
phase separation happens between states with the same 
magnetic ordering but different density. Densities ni and 
n2 are determined as above. Note, that while at n = 6 
DS is intrinsically unstable towards phase separation, at 
this density phase separation between DS and (tt, 0) is 
favored. 



